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Abstract 

Based on the compilation of the available e + e~ data Q, we present a non-perturbative 
estimation of the Adler function derived from the electromagnetic current correlator, and 
compare it with theoretical predictions from perturbative QCD (pQCD). The comparison 
is presented for the Euclidean region where pQCD is supposed to work best. We empha- 
size that such a comparison only makes sense if one takes into account the exact mass 
dependence of the perturbative predictions, which are available for the leading and next 
to leading (two-loop) order. In order to have the correct physical mass dependence in 
the evolution of the strong coupling as well, we utilize the MOM scheme /3-function to 
two-loops calculated recently Three-loop effects, which are available as series expan- 
sions for lowQ and highQ momentum transfer, are included by using Pade improvement. 
We discuss possible constraints on non-perturbative effects as suggested by the operator 
product expansion. 
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1. Introduction 

There are different ways of comparing QCD theoretical predictions with available time-like 
experimental data for the e + e~ total cross-section. Among the proposed methods are the 
global duality sum rules with different smearing functions [§, || 0, § as well as the local duality 
sum rules ||, [H| . The Borel sum rules have already been used before to analyze the low-energy 



e + e data [11] and gave the reasonable value of the gluon condensate < G > and a rather 



low value of the parameter A-^, which turns out to be in contradiction with the current world 



average value of a s (M z ) = 0.119 ± 0.002 |L2 



In view of the improved understanding of the influence of the perturbative QCD corrections 
(both massless [0, |14| and massive ones || f|, [15|]) and the appearance of the compilation of 
the available e + e~ data |l]| it is worthwhile to study the problem of the comparison of the 
experimental data with the theoretical predictions at the new level. 

In this our work we will follow the idea originally proposed in Ref . [[16] and implemented in part 



in Ref.|T7| and will use for the comparison with theory the experimental data for the Adler 



function, defined through the following dispersive relation 

where Q 2 = —q 2 is the squared Euclidean momentum transfer, s is the center of mass energy 
squared for hadron production in e + e~-annihilation, 

D/ v a tot (e + e~ -> 7* -> hadrons) 
v ' cr(e+e- -> 7* -> fM+fi') 

the hadronic to leptonic cross section ratio with cr(e + e~ — > 7* — ► = ^fj-, ct = e 2 /4ir 

being the QED fine structure constant. 

The attractive features of this procedure were advocated from different points of view in Ref. [18 
and Ref. [0. 



2. Calculations of the Adler function 

In this note we reconsider the question of the reliability of pQCD predictions at lower energies 
(down to about 1 GeV) by calculating the Adler function associated with the hadronic part of 
the electromagnetic current 

2 1 - 

q i 6 6 

By Ui and dj we denote the upper and lower components of the quark weak-iso doublets, re- 
spectively (z is the family index). 

is the conserved flavor diagonal vector current of the quark q with charge Q q and mass m q . 
We first define the photon vacuum polarization amplitude 11^ (q 2 ) by 

KM = if*™** < o\tj;(x) T> (0) [0 > 

= - (V 9nv - q^qv) K (q 2 ) ■ (3) 



2 



By U' v we denote the corresponding amplitude for the vector current J^y- 
The Adler function is defined as the derivative 



2 , dW^s) 



D(-s) = -(12tO a . (4) 



We may write 



where 



D{-s) = — s-^-Aa had (s) (5) 
a as 

Aa had (g 2 ) = Re f°° ds . R ^ , . (6) 

37T J Ami s \ s - Q l£ ) 

is the hadronic contribution to the shift of the fine structure constant at scale q 2 . A careful 
estimate of this function using the existing experimental data was given some time ago in 
Ref. Experimental data for R(s) may be used up to E cnt = 40 GeV, for higher energies 
7Z mixing would complicate the analysis. The high energy tail is calculated by using pQCD. 
Up to three-loops results are available for massive quarks |2(J. The four-loop corrections are 



known in the approximation of massless quarks [14 



In the present note we are going to evaluate the Adler function by the same procedure. It is 
important that we work in the Euclidean region Q 2 = —q 2 > where pQCD is supposed to 
work best. 

Before presenting the non-perturbative result, which follows from the estimate of the dispersion 
integral, we briefly discuss the analytic results which were obtained from perturbative QCD 
and the operator product expansion (OPE). We write 

D(Q 2 ) = D(°)(Q 2 ) + D {1) (Q 2 ) + D^ 2 \Q 2 ) + ■■■ + D NP (Q 2 ) (7) 

where D^ n \Q 2 ) is the n-th order pQCD contribution proportional to (a s /n) n and D NP (Q 2 ) is 
the non-perturbative part which will be specified below. 

We first consider the pQCD contributions. For the representation of the explicit expressions 
we use the variables 



2/ _ , VT^y-l 



imi/s , f 



where £ is taking values 0<£<lfors<0. 

The MS subtracted transverse photon vacuum polarization function 11^ (q 2 ), which defines the 
full photon propagator as 

i 1 



'q 2 l + e 2 W 1 (q 2 
at the one-loop level is given by 



with fi the MS scale parameter, N c f the color factor and 

y 



The corresponding one-loop contribution to Act is 

Aa(q 2 )=e 2 (ifyO) - Tl'^)) 

while for the one-loop Adler function we obtain 

D^(Q')=J2Q}N cf H^ (9) 
/ 

with 

H = (12tt 2 ) tl' v , tl' v = -s dU f v /ds , (10) 
in terms of the vector current amplitude W v . Explicitly we have 

y dy 2 4 q 

Asymptotically, we find 

H (0) ^ I ^ " ^ + f 5 + ' ' ' « ™J 

and this behavior determines the quark parton model (QPM) (leading order QCD) property of 
the Adler function: heavy quarks (mj Q 2 ) decouple like Q 2 /rrvj while light modes (m 2 <C Q 2 ) 
contribute Q 2 N c j to D^°\ Using the two-loop QCD vacuum polarization functions given in 
the appendix we obtain for the two-loop contribution to the Adler function 



D^\Q 2 ) = ^±Y,Q 2 fN cf H^ (11) 



7T / 



where = (12tt 2 ) tl v 2) (-Q 2 , m}) . 

The full massive three-loop vector current correlator n^ 3 -* is not yet known, however, a sufficient 
number of terms were calculated both for the small mass expansion @| as well as for the heavy 
mass expansion ||. By taking the derivative (^) we obtain the contribution to the Adler 
function 

rt 2 \Q 2 ) = (^p-Jj2Q}N cf H^ (12) 
where = (I2ir 2 ) fl v 3) (-Q 2 , m 2 f ) . 

The convergence of the series expansions break down at Euclidean Q 2 = 4m 2 , right in the 
region where mass effects are most important, and we use a conformal mapping together with 
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Pade improvement [^TJ in order to obtain results useful for the present analysis. For details we 
refer to a brief discussion in the appendix. 

A way to parametrize non-perturbative (NP) effects away from resonances, as applicable here, 
is prescribed by the OPE of the current correlator (|3D, valid for large enough Q 2 . This part 
is due to non- vanishing gluon and light quark condensates 0j. They yield the leading power 
corrections 



D NP (Q 2 



q=u,d,s 



8tt z 



'-GG > 



+ 



+ 




(13) 



+ 



where l m = ln(Q 2 //z 2 ). < ^-GG > and < m q qq > are the scale-invariantly defined condensates. 
The terms beyond leading order in a s have been calculated from the results which were obtained 
in Refs. [^3], |[23|| , p4]| and [fZIJ. The terms beyond the leading order in a s are unimportant 
because experimental data only yield rough constraints on the condensates and we include them 
for completeness only. Sum rule estimates of the condensates yield |7|, f26| 



< —GG >~ (0.336 ■ ■ ■ 0.442 GeV) 4 and < m q qq >'. 

71 



(0.086 • 
-(0.192 • 



•0.111 GeV) 4 for q 
•0.245 GeV) 4 for q 



u, d 



with rather large uncertainties. The values assumed for the quark condensates are based on 
the recent estimate|27] < qq >~ — (235 ± 27 MeV) 3 and the quark masses m u 
(m u + m d )/2 ~ 7.2 ± 1.1 MeV and m s 
be evaluated numerically at /i 2 = Q 2 . 



175 ± 25 MeV pi, all at scale 1 GeV. Eq. (|T3|) will 



3. Numerical results and conclusions 



For numerical studies we used the following pole quark masses |T2 

m u ~ ma ~ m s ~ ; m c = 1.55GeV ; nib = 4.70GeV ; m t = 173.80GeV . 
For the strong interaction coupling we take «^^g = 0.120 ± 0.003 at the scale M^=91.19 



GeV 



of the amplitudes as 



] . Up to two-loop we use the exact mass dependence (see e.g. 
well as of the exact mass dependent two-loop a s (Q 2 ) in the gauge invariant background-field 
MOM renormalization scheme, presented recently in Ref. 0. In the transition from the MS to 
the MOM scheme we adapt the rescaling procedure described in 0, such that for large [i 

a s ((x /i) 2 ) = a s (/j 2 ) + + 0(a 3 s ) . 

This means that xq is chosen such that the couplings coincide to leading and next-to-leading 
order at asymptotically large scales. Numerically we find xq — 2.0144. Due to this normaliza- 
tion by rescaling the coefficients of the Adler-function remain the same in both schemes up to 
three-loops. 
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The numerical evaluation of the dispersion integral (|1|) is based on the compilation of data 
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Figures la and lb: We show the "experimental" Adler function together with pQCD pre- 
dictions and power corrections (NP) in the region below about 10 GeV (a) and separately 
for the low energy tail (b). The negative sign chosen for Q should remind the reader that 
we are considering space-like momentum transfer. The shaded area represents the ±lcx 
band obtained from the data. A similar error band is shown for the uncertainties in the 
power corrections. 

discussed in [ED]. Here we will not use all the data up to 40 GeV, however, still take a con- 
servative attitude, and replace data by pQCD results only where it is obviously safe to do so. 
Accordingly, in the regions from 4.5 GeV to Mx and above 12 GeV we use pQCD, including the 
massive results available up to three-loops [f| and the four-loop contribution in the massless 
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approximation [13]. The "experimental" behavior of the D-function is displayed in Figs. 1(a) 
and 1(b), together with the theoretical predictions. The data and the error analysis utilized 
in the evaluation of (|I|) are described in Ref. [|I|]. Both statistical and systematic errors of the 
data were taken into account. Note that the meaning of the la error band differs from the 
usual "one standard deviation" since its points at different Q are highly correlated. The origin 
of the uncertainty is shown in Table 1. 

Table 1: Origin of uncertainties for the "experimental" Adler function at Q=2.5 GeV and Mz- 
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At Mz D(Q 2 ) is completely dominated be the perturbative high energy tail, the uncertainty is 
the one of pQCD. As Q 2 decreases the dominating region moves to lower energies. At 2.5 GeV 
contributions to D(Q 2 ) are spread over all energies below Mr and uncertainties from the region 
below Mj/^ are dominating. Errors at 2.5 GeV are about a factor 20 larger and dominated by 
systematic errors of the e + e~ data. 

The agreement between pQCD and the "data" looks rather good, but only after inclusion of the 
3-loop contribution. If we had not include the 3-loop term it would appear that there is some 
room for the non-perturbative power corrections (|13|). After adding the 3-loop result there is 
no region left where the power corrections, with the above given values for the parameters lead 
to an improvement. So, the comparison of the perturbative predictions with the "data" looks 
quite impressive after inclusion of the 3-loop massive contribution. The inclusion of the top 
(nf=6) does not alter the results at the energies shown. We can not see a region of intermediate 
Q 2 where adding the non-perturbative power corrections (NP) does improve the agreement of 
3-loop pQCD predictions with the presented "data". It also does not look possible to estimate 
in a convincing way the leading coefficient < ^-GG >. Due to the power behavior, these kind 
of corrections drop very fast at larger energies, above about 2 GeV, while they grow very fast 
at lower energies, below about 0.75 GeV, which signals the breakdown of the applicability of 
the expansion. In fact, including the l/Q 6 terms would lead to a large contribution of opposite 
sign at low Q 2 . The main problem is that these corrections become substantial only at energies 
where the perturbation expansion can no longer be trusted. 

It should be noted that in the Euclidean region there are no "threshold steps" as they are 
suggested if the MS scheme is applied. In the MS scheme effective theories with different 
number of light flavors must be matched at the different scales set by the quark masses BIJ in 



order to restore decoupling of heavy flavors. The physical behavior is illustrated in Fig. 2. 
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Figure 2: The flavor dependence of the Adler function in pQCD. For each nf = 3 • • • 6 two 
curves are shown: one for the QPM (lower of the two) and the other for pQCD at three-loops 
(upper of the two). The cases nj =5 and nj = 6 are graphically indistinguishable at the 
energies shown. The thickness of the n/ = 5 0(a 2 ) "line" illustrates the small theoretical 

(5) 

uncertainty obtained from the variation a = 0.120±0.003 at the scale M^=91.19 GeV. 



In particular at low momenta the MS prescription does not provide a good description of 
the physical behavior and hence its applicability in the Euclidean region is worth a critical 
inspection. In this spirit we also would expect an improvement of the QCD description of 
deep inelastic scattering in cases where different effective numbers of flavors n/ have to be 
considered 



Note that the 3-loop term is quite large and higher orders still could give sizable contributions. 
The massless 4-loop contribution may be estimated as follows: in massless pQCD we have 
the familiar result [HJ] R(s) = x (1 + a + cia 2 + C2a 3 + • • •) with a = a s (s)/ii, c\ 



1.9857 - 0.115372/ and c 2 = -6.6368 - 1.2002/1/ - 0.0052^ - 1.2395 (E<3/)7( 3 £<2/)- In 
this approximation R{s) is a constant if we discard the running of a s and we can evaluate (jl]) 
obtaining D(Q 2 ) ~ R{Q 2 ) at large Q 2 . Numerically the 4-loop term proportional to c 2 amounts 
to —0.06% at 100 GeV and increases to about —0.5% at 2.5 GeV. At the low energy this estimate 
is of course not reliable because we do not expect the approximation to work. Nevertheless in 
may indicate the possible size of the missing higher order corrections. Due to the increase of 
the running strong coupling the higher order predictions become unreliable at lower energies. 
Once we observe that the functional dependence deviates from that of the data, which happens 
below about 2.5 GeV, pQCD is not reliable any longer. The large uncertainties in the data (see 
Table 1) on the one hand and the beginning failure of pQCD on the other hand, unfortunately 
do not allow to draw any conclusion about non-perturbative physics. Calculations of missing 
higher order terms including the mass effects certainly would help a lot to draw more reliable 
conclusions about the borderline of validity of pQCD. 
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Appendix 

Here we present the explicit formulae for the two-loop vacuum polarization amplitudes of the 
flavor diagonal vector current. We first introduce the abbreviations 

/ = -~ln£, </ = ln(l-£), h = ln(l + £) 

where y and £ are given by flD and 

ALi 3 = 2Li 3 (0 - Li 3 (£ 2 ) , ALi 2 = Li 2 (0 - Li 2 (^ 2 ) 
and define the auxiliary functions 

XX = ALi 3 + ^/ALi 2 + ^/ 2 (2/i + £), 
YY = ^(ALi 2 + 2/(2/i + s) + 3/ 2 ) . 
The poly-logarithms IA n (z) are defined by the integrals 

Un{z) ~ (^2)! Jo x dX 

and ( n = ((n) = X^il/ 77 -* is the Riemann zeta-f unction, with values ( 2 = 7I " 2 /6, C3 = 
1.202057..., (4 = vr 4 /90 etc.. 

In terms of these functions we find the following compact expression for the two-loop contri- 
bution to the Adler function related to the QCD vector neutral current vacuum polarization 
amplitude:! 2 ] 

127r 2 n^ 2) = -2 2/ 2 XX+ v /T^((l + 2/ )YY-14/3 2 //) 
+ (3y 2 - 4 - 4/(1 - y)) f 
+ l + ll/3 2 / + 2y 2 C(3). 

The following exercise about handling asymptotic expansions is intended to illustrate at the 
two-loop level the procedure to be applied at the three-loop level where an exact massive result 
is not yet known. 



2 it is given by 

I^Ul (Re< c ( s )/ s -ReITr( S )) 

in terms of the amplitudes which we can find in Sec. 4.5 of Ref. (see also |p3||). There the rules for the 
analytic continuation to the time-like region can also be found. 
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The expansion for large Q 2 (y = —4m 2 /Q 2 — > 0) reads ( L = \n(Q 2 /m 2 ) ) 



i27r 2 rr'f 
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+ 
+ 
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and for low Q 2 we have 
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Both series expansions diverge at the boundary of the circle of convergence Q 2 = 4m 2 . In fact 
this is precisely in the region where mass effects are of the order of unity in the Euclidean 
region. We therefore apply a conformal mapping 
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Figure 3: The two-loop amplitude H(l) = (127r 2 )ily 2 ^: shown is the exact result (dotted 
line), the series expansions for low Q 2 (8 terms) (4m 2 > Q 2 ) and for large Q 2 (7 terms) 
(Q 2 > Am 2 ) (dashed line) and the Pade improvement (full line) of the low Q 2 series working 
excellent up to (16m 2 > Q 2 ). For higher Q 2 the large Q 2 -expansion works perfectly. The 
dotted line is not seen because of the perfect agreement between the exact result and the 
Pade improved approximation. 
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from the complex negative q 2 half-plane to the interior of the unit circle \oj\ < 1 together with 
Pade resummation as proposed in [^TJ. The Pade approximant provides a good estimation 
to much higher values of 1/y up to about 1/y ~ 4. This is displayed in Fig. 3, which shows 
that utilizing the Pade improvement allows us to obtain reliable results also in the relevant 
Euclidean "threshold region", around y = 1. 

The three-loop vector current amplitude II ( 3 ) is available in the form of a small Q 2 expansion, 
Eq. (47) of Ref. [§[, and as a large Q 2 expansion, Eqs. (7-10) of Ref. |§. By the definition (|) 
the contribution to the Adler function is obtained by differentiation 

tl'^(y) = -ydU'^(y)/dy . 

While the low Q 2 expansion is given in the on-shell scheme, where it is a simple power series in 
1/y, the high Q 2 expansion is given in the MS scheme as an expansion in y = 4m 2 /Q 2 where 
m = m (/i 2 ) is the MS running mass at scale /i and the MS amplitude IT'® includes powers of 
the logarithms l qm = ln(— q 2 /m 2 ) and l qfl = ln(— q 2 / fi 2 ) in additions to powers of y. 

Since, for simplicity, we work with pole masses, we have to reparametrize the MS amplitude. 
The relationship between the MS mass fh and the pole mass m is given by p, [34], [H| 

— ( 1+ -(s) + -(s) V -)' < 14 > 

with 

ci = C F (4 + 3L), (15) 

c 2 = C F C A (HH _ 8C2 _ 4/3 (l) + + ^L 2 ) 

/71 26 \ 

- C F T F n f (- + 8C2 + -L + 2L 2 j 

/ 1 21 27 Q \ 

+ c 2 F (— + 30(2 + 8/3(1) + 2 L J ~ 12CfTf - 2 & ■ ( 16 ) 

where / 3 (1) = |Cs - 6C2 In 2, L = ln(/i 2 /m 2 ), C A = 3, C F = 4/3, T F = 1/2 and n f = 6. 
Asymptotically, for large Q 2 we find (l qil = ln(— q 2 /fi 2 )) 

H^ioo) = ^{(-3C 2 + (123 - 88 C 3 - 22 Z W )C A C F 

+ (-44 + 32 C 3 + 8 V) C F T F m + (-44 + 32 Cs + 8 l qil ) C F T F )} 

which coincides with the MS expression for the a 2 coefficient in the D-function, calculated in 
Ref. At the scale /1 2 = Q 2 is of order unity: H^ 2 \oo) ~ 1.40923. 

Due to the logs ln(Q 2 /m 2 ), still present after setting /z 2 = Q 2 in the three-loop amplitudes, 
we perform the Pade improvement for the non-log and the log terms separately. The result is 
depicted in Fig. 4. 
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Note that the Pade approximant is an analytic continuation of the small Q 2 expansion with 8 
terms. From a practical point of view it provides the "exact" mass dependence of the three-loop 
amplitude. 
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Figure 4: The three-loop amplitude = (127r 2 )iTy 3 ' ) : shown are the series expansions 
for low Q 2 (8 terms) (4m 2 > Q 2 ) and for large Q 2 (7 terms) (Q 2 > 4m 2 ) (dashed line) and 
the Pade improvement (full line) of the low Q 2 series working excellent up to (16m 2 > Q 2 ). 
For higher Q 2 the large Q 2 -expansion works perfectly. 
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